SOME ERROR ESTIMATES FOR THE FINITE VOLUME 
ELEMENT METHOD FOR A PARABOLIC PROBLEM 



P. CHATZIPANTELIDIS , R.D. LAZAROV, AND V. THOMEE 

j^JQf Abstract. We study spatially semidiscrete and fully discrete finite volume el- 

r ^ _ ement methods for the homogeneous heat equation with homogeneous Dirich- 

let boundary conditions and derive error estimates for smooth and nonsmooth 
initial data. We show that the results of our earlier work \5\ for the lumped 
mass method carry over to the present situation. In particular, in order for 
error estimates for initial data only in L2 to be of optimal second order for 
positive time, a special condition is required, which is satisfied for symmetric 
triangulations. Without any such condition, only first order convergence can 
be shown, which is illustrated by a counterexample. Improvements hold for 
triangulations that are almost symmetric and piecewise almost symmetric. 



1. Introduction 
We consider the model initial-boundary value problem 

(1.1) u t — Au = 0, in il, u = 0, on dfl, for t > 0, with u(0) = v, in £1, 

where fHs a bounded convex polygonal domain in K 2 . We restrict ourselves to the 
homogeneous heat equation, thus without a forcing term, so that the initial values 
v are the only data of the problem. This problem has a unique solution u(t), under 
appropriate assumptions on v, and this solution is smooth for t > 0, even if v is 
not. 

To express the smoothness properties of the solution of let, for q > 0, 

H q C £2(^2) be the Hilbert space defined by the norm 

(1.2) \w\ q — {^2/ ^( w ^j)j ' where (w,<p) = / wipdx, 



i=i 



and where {Xj}'^. 1 , {4>j}™^ l are the eigenvalues, in increasing order, and orthonor- 
mal eigenfunctions of —A in ft, with homogeneous Dirichlet boundary conditions 
on <9f2. Thus \w\o = \\w\\ = (w, w) 1 ^ 2 is the norm in L 2 = L 2 (Q), \w\\ = ||Vw|| 
the norm in Hq = Hq(Q) and |to|2 = is equivalent to the norm in H 2 (il) 

when w — on dfl. Eigenfunction expansion and Parseval's relation shows for the 
solution u(t) = E(t)v of (|1.1|> the stability and smoothing estimate 

(1.3) \E(t)v\ p <Ct-^ p - q ^ 2 \v\ q , for0<g<p, and t > 0. 

In fact, since the smallest eigenvalue is positive, a factor of e~ ct , with c > 0, may 
be included in the right hand side, and this holds for all our stability, smoothing 
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and error estimates throughout our paper. Since our interest here is in small time 
we shall not keep track of this decay for large time below. We shall also use the 
norm = J2\-y\<k su Pzen \^I W ( X )\ m C k = ^ fc (^)i wrtn C = the space of 

continuous functions on Here for 7 = (71,72), T>1 = (d/dxi)' 11 (d/dx 2 ) 12 and 
|7| = 71 + 72- 

We first recall some facts about the spatially semidiscrete standard Galerkin 
finite element method for (|1.1[) in the space of piecewise linear functions 

S h = {x^C: x\t linear, Vr e %; x\dQ = 0}, 

where {Th} is a family of regular triangulations Th — {t} of f2, with h denoting the 
maximum diameter of the triangles r £ Th- This method defines an approximation 
Uh(i) G 5/j of u(t), for i > 0, from 

(1.4) (u M , x) + (Vufc, Vx) =0, Vx e S h) for t > 0, with u h (0) = v h) 

where Vh £ Sh is an approximation of v. It is well-known that we have the smooth 
data error estimate, valid uniformly down to t — 0, see e.g. |15j , 

(1.5) \\u h (t) - u(t)\\ < Ch 2 \v\ 2 , if \\v h -v\\ <Ch 2 \v\ 2 , fort>0. 

We also have a nonsmooth data error estimate, for v only assumed to be in L 2 , 
which is of optimal order 0(h 2 ) for t bounded away from zero, but deteriorates as 
t -)• 0, 

(1.6) \\u h (t) - u(t)\\ < CTi^lMI, if v h = P h v, for * > 0, 

where Ph denotes the orthogonal L 2 — projection onto Sh- Note that the choice of 
discrete initial data is not as general in this case as in (|1.5p . We emphasize that 
the triangulations Th are assumed to be independent of t, and thus the use of finer 
Th for t small is not considered here. 

We note that a possible choice in (jl.5|) is Vh — PhV, and hence, by interpolation, 
we have the intermediate result between (|1.5p and (ll.6[) . 

(1.7) \\u h (t) -u(t)\\ < Chh'^lvlu i£v h =P h v, fort>0. 

Recently, in [5], we showed results similar to (|1.5p - (|1.7[) for the lumped mass 
finite element method, which may be defined by replacing the L 2 — inner product 
in the first term in (|1.4[) by the quadrature approximation (uh.t,x)h, where, with 
Ih '■ C — > Sh being the interpolant defined by Ihv(z) — v{z) for any vertex z of Th, 

{X,i>)h = / Ih{xi>)dx, \/x,i> € Sh. 
Jn 

Improving earlier results, we demonstrated that (| 1 . 5|) remains valid for the lumped 
mass method, but that (jl.6p requires restrictive conditions on {Th}, caused by the 
use of quadrature in (jl.4p . and satisfied, in particular, for symmetric triangulations. 
We remark that the choice of discrete initial data in the analogue of (|1.7p was 
incorrectly stated in [5] , see Section [3] below. 

In the present paper our purpose is to carry over the analysis in [S] to the 
finite volume element method for problem ([l.ip . This method is based on a local 
conservation property associated with the differential equation. Namely, integrating 
(jl.ip over any region V C and using Green's formula, we obtain 

(1.8) / u t dx- Vu-nda = 0, for t > 0, 

JV JdV 
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Figure 1. Left: A union of triangles that have a common ver- 
tex z; the dotted line shows the boundary of the corresponding 
control volume V z . Right: A triangle r partitioned into the three 
subregions t 2 . 



where n denotes the unit exterior normal vector to dV. The semidiscrete finite vol- 
ume element approximation Uh(t) £ Sh, will satisfy (|1.8[) for V in a finite collection 
of subregions of f2 called control volumes, the number of which will be equal to the 
dimension of the finite element space Sh- These control volumes are constructed 
in the following way. Let z T be the barycenter of r £ Th- We connect z T by line 
segments to the midpoints of the edges of r, thus partitioning r into three quadri- 
laterals t z , z £ Z/j(t), where Z%{j) are the vertices of r. Then with each vertex 
z £ Zh = yj T £j- h Zh{T) we associate a control volume V z , which consists of the union 
of the subregions r z , sharing the vertex z (see Figured! left). We denote the set of 
interior vertices of Zh by Z^. The semidiscrete finite volume element method for 
(|l.lj) is then to find Uhif] £ Sh such that 

(1.9) / Uh^tdx— I Vuh ■ nda = 0, Vz £ ZjJ, for t > 0, with Uh(0) = Vh, 
Jv z JdVz 

where Vh G Sh is an approximation of v. 

This problem may also be expressed in a weak form. For this purpose we intro- 
duce the finite-dimensional space of piecewise constant functions 

Y h = {v G L 2 : r)\ Vz = constant, £ ?y|vi =0, Vz e Z h \ Z^}. 

We now multiply (|1.9|) by t?(z) for an arbitrary 77 £ Y/j, and sum over z £ Z2 to 
obtain the Petrov-Galerkin formulation 

(1.10) (Uh,t,v) + ah(uh,v) = 0) V?7 G Y^, for i > 0, with S A (0) = v h , 
where the bilinear form a^(-, •) : Sh X Yh — > R is defined by 

(1.11) a h (x,r]) = - r)(z) Vx-nda, V% £ S/,, 77 £ F^. 

Obviously, we can define a^(-, •) also for % replaced by u> £ H 2 , and using Green's 
formula we then easily see that 



Oh (to, ?y) = -(Aw, 17), Vw £ iJ 2 , 77 £ Y h - 



4 



ERROR ESTIMATES FOR THE FVEM 



We shall now rewrite the Petrov-Galerkin method (|1.10j) as a Galerkin method 
in Sh- For this purpose, we introduce the interpolation operator J/j : C t— > Yh by 

JhU = ^ u ( z )*2' 

where \& 2 is the characteristic function of the control volume V z . It is known that 
Jh is selfadjoint and positive definite, see [6], and hence the following defines an 
inner product (•, •) on Sj,, 

(1.12) (x,^) = (x,AV0, Vx,^e5 h . 

Also, the corresponding discrete norm is equivalent to the L2— norm, uniformly in 
/i, i.e., with C > c > 0, 

c||x|| < Hlxlll < C||x||, Vx 6 S h , where ||| x ||| = (x,x} 1/2 , 

see [Sj. Further, in [2 , it is shown that 

ah{x,Jhip) = (Vx,VV>), Vx,^eS/„ 

and therefore, a/i(-, •) is symmetric and a%[x, JhX) — ll^xll 2 ; f° r X € Sti- 
With this notation, (| 1 . 10[) may equivalently be written in Galerkin form as 

(1.13) («MiX) + (V5fc,Vx)=0, VxeSft, fort>0, with u h (0) = v h . 

Our aim is thus to show analogues of (|1.5|) - (|1.7|) for the solution of (| 1 . 1 3|) . with 
the appropriate choices of Vh, i-e., 

(1.14) \\u h {t)-u{t)\\ <Ch 2 t- 1+q l 2 \v\ q , fort>0, g-0,1,2. 

This will be done below for q = 2, and in the case q = 1 under the additional 
assumption that {%,} is quasiuniform. However, for q — 0, as in [5], we are only able 
to show ([1.14)1 under an additional hypothesis, expressed in terms of the quadrature 
error operator Qh '■ Sh —> Sh, defined by 

(1.15) (VQ /l V,Vx) = £/ [ (V,X), Vx^G^, 
where £/ l (-, •) is the quadrature error defined here by 

(1.16) e h (f, X ) - (/, JhX) - (f,x), V/ G L 2 , x e S fc , 
and requiring 

(1.17) HQfcVll < C7* a M, V^eSh. 

We will show that this assumption is satisfied for symmetric triangulations Th- 
Symmetry of Th, however, is a severe restriction which can only hold for special 
shapes of Q,. For this reason we will also consider less restrictive families {Th}- 
We will demonstrate that (| 1 . 1T[) holds for almost symmetric families (discussed 
in Section [4]), with the addition of a logarithmic factor; we also show that this 
logarithmic factor is not needed in one space dimension. Further, for piecewise 
almost symmetric families of triangulations, see Section 21 the inequality (|1.17|) 
holds with an 0{h 3 / 2 ) bound. 

We then give two examples of nonsymmetric triangulations such that ([1. 141) 
does not hold for q = 0. In the first example we construct {Th} such that the 
convergence factor is at most of order 0(h) for t > 0, and in the second example, 
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with nonsymmetry only along a line, of order 0(h 3 ^ 2 ). Without any additional 
condition on Th we are only able to show the nonoptimal order error estimate 

\\u h {t) -u{t)\\ < Cht- 1/2 \\v\\, i{v h = P h v, fori>0. 

We remark that in p2] , in the more general case of a parabolic integro-differential 
equation, the nonsmooth data error estimate (|1.14p . for q = 0, with an extra factor 
| log/i|, was stated, for any quasiuniform family {Th}- Unfortunately, this result is 
in contradiction to our above counterexamples, and its proof incorrect. 

We also discuss optimal order O(h) error estimates for the gradient of Uh — u, 
under various assumptions on the smoothness of v and choices of Vh- Further, in 
a separate section, we consider briefly the extension of our results for the spatially 
semidiscrete problem to the fully discrete backward Euler and Crank-Nicolson finite 
volume methods. 

As for the lumped mass method in [5] , our analysis yields improvements of earlier 
results, in 3 >; where it was shown that, for smooth initial data and Vh = RhV, 

\\u h (t) - u(t)\\ <Ch 2 \v\ 3 , fort>0, 

and 

||V(SfcO) - u(t))\\ < Che^lv^+e, for <> 0, e > small. 

As in the case of the lumped mass method in [5], these improvements are made 
possible by combining, the error estimates (|1.5p - (|1.7l) for the standard Galerkin 
finite element method with bounds for the difference 8 — Uh — Uh, which, by (|1.13p 
and (|1.4[) . satisfies 

(1.18) (SuX) + (V5,V X ) = -£h(u h ,t,x), Vxe5 ft , for t > 0. 

In the final section we sketch the extension of the theory developed above to 
more general parabolic equations, considering the initial-boundary value problem 

(1.19) u t + Au = 0, in fi, u = 0, on dtt, for t > 0, with u(0) = v, in Q, 

where Au — —V ■ (aVu) + /3m, with a a smooth symmetric, positive definite 2x2 
matrix function on £1 and /3 a non-negative smooth function. 

Here, let Uh(t) S Sh, denote the standard Galerkin finite element approximation 
of u(t), defined by 

(1.20) (u h ,t,x)+a(uh,x) = 0, VxeS 1 ^ for t > 0, with u h (0) = v h , 
where Vh G Sh is an approximation of v and 

(1.21) a(u>, ip) — (aVto, V<p) + {(3w, if), for w,ip £ Hq. 

In a straight -forward way the estimates (|1.5|) - (|1.7|) extend to the solution of (|1.20[) . 

The natural generalization of the finite volume method (|1.10[) would now be to 
find Uh(t) G Sh such that 

(1.22) (u h ,t,x) + a h (u h , JhX) = 0, Vx £ for t > 0, with u h (0) = v h , 
where, instead of (|1.11|) . one uses the bilinear defined by 

(1.23) a h (fl>,v) = 7 ?( z )(- / (a^)-nda+ / Ptj)dx\, € Sh, f) S Yh- 
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It is known that, in general, the bilinear form ah(tp, JhX)i is nonsymmetric on Sh but 
it is not far from being symmetric, or a/i(x, Jh4>) — ah(*P, Jhx)\ < C^ll^xll IIW>||, 
cf. [BJ. Also, if a and j3 are constants over each r G 7ft, then, see, e.g. [HIH], 

(1.24) a h {il>, JhX) = («V^, V X ) + (P1>, JhX), W>, X € S h , 

and thus a>h(ip,JhX) 1S symmetric, since as we shall show (flip, JhX) = (PXi Jh^)- 
Therefore, since symmetry is important in our analysis, we introduce the modified 
bilinear form 

(1.25) a h (ip,v) = Y\ V(z)(~ I (aV^)-nda+ [ ffldx), W> G S h , r, G Y h , 

where, for z G r, r G 7h, a(z) = a(z T ) and j3{z) — /3(z T ), with z T the barycenter 
of r. This choice of a^(-,-) leads to the finite volume element method, to find 
Uh{t) G Sh such that 

(1.26) (uh,t,x) + ah(u h , Jhx) = 0, V% G SVi, for i > 0, with 5^(0) = V/j, 

and for this the desired analogues of the estimates (jl.l4j) are established in Theo- 
rems EEO 

The following is an outline of the paper. In Section [51 we introduce notation 
and give some preliminary material needed for the analysis of the finite volume 
element method. Further, we derive smooth and nonsmooth initial data estimates 
for the gradient of the error in the standard Galerkin method. In Section [3] we 
derive the error estimates (jl.141) discussed above, under the different assumptions 
on smoothness of data and the triangulations {Th}- In Section [4] we show that 
assumption (|1.17j) is valid for symmetric meshes, and discuss the corresponding 
properties for almost symmetric and piecewise almost symmetric meshes. In Section 
[S] we present two nonsymmetric triangulations in two space dimensions for which 
optimal order L2 _ convergence for nonsmooth data does not hold. In Section [6J 
we consider briefly the application to the fully discrete backward Euler and Crank- 
Nicolson finite volume methods. Finally, Section[7]contains the extension of Section 
[3] to more general parabolic equations. 



2. Preliminaries 

In this section we show a smoothing property for the finite volume element 
method, and discuss the quadrature associated with this method. We also derive 
some estimates for the gradient of the error in the standard Galerkin finite element 
method which will be needed later. 

We first recall that for the standard Galerkin method, one may introduce the 
discrete Laplacian : Sh — >• Sh by 

-(A^, x ) = (v^,Vx), WsxeSfc, 

and write the problem (|1.4|) as 

(2.1) u h ,t~&hUh = 0, for t > 0, with u h (0) =v h - 

Letting {A^}^ , {(fij}^ , where Nh = dim Sh, denote the eigenvalues, in increasing 
order, and the corresponding eigenfunctions of —Ah, orthonormal with respect to 
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(•,•), we have for the solution operator Eh(t) — e Aht of ()2.1|) . by eigenfunction 
expansion, 

N h 

u h (t) = E h {t)v h = XV X '*(«A>$)$. for t > 0. 

The following smoothing property analogous to (|1.3I) holds for v h £ Sh and t > 
(2.2) ||V^E h (tK|| <C^-^/ 2 ||V^ fe ||, £>0, p,<? = 0,l, 2£ + p>q, 
with A = d/dt. 

Turning to the finite volume method (|1.13[) , we now introduce the discrete Lapla- 
cian Ah : Sh Sh, corresponding to the inner product (•, •) in (|1.12[) . by 

(2-3) -(A h V,x) = (VV,Vx), V^, X eS k . 

The finite volume method (| 1 . 1 3[) can then be written in operator form as 

(2.4) u h ,t-AhUh = 0, fori>0, with u h (0) = v h . 
For the solution operator Eh(t) — e Aht of (|2.4j) we have 

N h 

(2.5) u h {t) = E h {t)v h = J2 e~ A ''(«fc, for t > 0, 

3=1 

where {X 1 -}^ and {^}^ are the eigenvalues, in increasing order, and the corre- 
sponding eigenfunctions, orthonormal with respect to (•, •), of the positive definite 
operator — A^. For Eh(t) the following analogue of (|2.2j) holds, cf. [5J Lemma 2.1]. 

Lemma 2.1. For Eh defined by (|2.5I) we have, for Vh G Sh and t > 

||V p £>f^(*HII <Ct-^ p - 9)/2 ||V^|U>0, p,g = 0,l, 2£ + p>g. 

Proof. Introducing the square root G/j = (— A/j) 1 / 2 : 5/j — » S/j, of — A^, we get 

\\Vvhf = ((-Ah)v h ,vh) =Y,~ Xt jM H i) 2 = III^IH 2 - 

i=i 

Since the norms ||| • ||| and || • || are equivalent on Sh we find for t > 

N h 

\\VPD e t E h (t)v h \\ 2 < C\\\GlD l t E h {t)v h \\\ 2 = C^(^) 2 ^-^- 2A ? t (A^)«( % , <^> 2 

3 = 1 

< Ct-^ e+ P-^\\\Glvh\\\ 2 < Ct-^ 2e+ P-^\\S/ q v h \\ 2 . □ 

The quadrature error functional Eh{-, •) defined by (I1.16P has an important role 
in our analysis below. For this reason we recall the following lemma, cf. [3J. 

Lemma 2.2. For the error functional £h, defined by (|1.16p . we have 

|£ fe (/,V)l<C^ +9 ||VP/||||V^||, VfEH\i>eS h) and p,? = 0,l. 

Proof. Since J (Jhtp — ip) dx = for ^ linear in r, for any tgTj, see 0, we have 
that Jhip — is orthogonal to 5^, the set of piecewise constants on 7^. Hence 

V>) = (/, ^ - = (/ - iV, ^ - 1*0, 

where is the orthogonal projection onto Sh- The lemma now easily follows since 
|[^^ - -0H < C^HV^II and 11^/ - /|| < □ 
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The following estimate holds for the quadrature error operator Qh in (|1.15jl . 
Lemma 2.3. Let Ah and Qh be the operators defined by (|2.3|) and (1 1 . 1 5[) . Then 

\\VQ hX \\ + h\\A h Q hX \\ < ChP +1 \\V p xl V X eS h)P = 0, 1. 

Proof. By (|1.15p and Lemma r2.2[ with ip = QhX an d Q = 1, it follows easily that 

l|VQ, lX || 2 = e h (x, QhX) < ChP +x \\W X \\ \\VQ hX \\, for p = 0, 1, 

which shows the desired estimate for ||V<3ftXll- Also, by the definition of Ah, 
Lemma |2~21 with q = shows, for p = 0, 1, 

\W&hQhX\\\ 2 = -(VQ hX ,VA h Qhx) = -e h ( X ,A h QhX) < ChP\\VP X \\ ||A h Q h x||. 
Since the norms ||| • ||| and || • || are equivalent on Sh, this implies the bound for 
remaining term || AhQhXW- ^ 

In addition to the orthogonal L2-projection Ph, our error analysis will use the 
Ritz projection Rh : Hq —> Sh, defined by 

(VR h w, V X ) = (Vw, V X ), V X e S h . 
It is well-known that Rh satisfies 

(2.6) \\R h w-w\\ +h\\V(R h w-w)\\ < Ch q \w\ q , for w € H q , q = 1, 2. 

We close with some estimates for the gradient of the error, slightly generalizing 
those of Theorem 2.1]. 

Theorem 2.1. Let u and Uh be the solutions of (|l.lj) and (|2.1j) . Then, for t > 0. 

(Ch\v\ 2 , if \\V(v h -v)\\ <Ch\v\ 2 , 

\\V(u h (t) - u(t))\\ < I Cht-^Mu if K - «|| < Ch\v\ u 
[Cht^Wvl zfv h = PhV. 

Proof. In [5J Theorem 2.1] this was shown with v% = RhV in the first two estimates, 
and thus it remains to bound V Eh(t)(vh — Rhv) ■ With d := Vfr—RhV we find easily, 
byLemmaO for smooth data, \\V E h (t)§(Q)\\ < |]Vi?(0)|| < Ch\v\ 2 , and for mildly 
nonsmooth data, \\VE h (t)d{0)\\ < Ct~ 1 ^\\'d(0)\\ < Ct^^h^. □ 

3. Smooth and nonsmooth initial data error estimates 

In this section we derive optimal order error estimates for the finite volume 
element method (|1.13j) . with initial data v in H 2 , H 1 and L 2 . For v 6 H 2 , the 
error estimate is the same as that for the standard Galerkin finite clement method, 
and this is also the case for v G H 1 . provided the family of finite element spaces is 
quasi-uniform. In the case v € L 2 , with discrete initial data Vh = PhV, in order to 
derive an optimal order estimate analogous to (|1.6p . we need to impose condition 
(|1.17p for the quadrature error operator Qh- In Section [4] we verify this condition 
for symmetric meshes. In the general case we are only able to show a non-optimal 
order 0(h) error bound in L 2l whereas for the gradient of the error an optimal 
order 0(h) bound still holds. 

The estimates and their proofs are analogous to those for the lumped mass 
method derived in [5], since the operators Eh, Ah and Qh, defined in Section [2l 
have properties similar to those of the corresponding operators for the lumped mass 
method. References to [5] will therefore be given in some of the proofs below. We 
begin with smooth initial data, v <G H 2 . 
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Theorem 3.1. Let u and Uh be the solutions of (jl.ljl and (|2.4j) . Then 
\\u h (t)-u(t)\\<Ch 2 \v\ 2 , if \\v h -v\\ <Ch 2 \v\ 2 , fort>0. 



Proof. Since, by (|1.5|) . the corresponding error bound holds for the solution Uh of 
the standard Galerkin method, it suffices to consider the difference 6 = Uh — Uh- 
Also, by the stability estimates of Lemma |2.1[ we may assume that u/j = RhV. By 
the definition (|1.15j) of Qh, S satisfies (jf . 18[) . and hence 

(3.1) S t - A h S = AhQhUh,t, for*>0, with 5(0) = 0, 
where Uh is the solution of fj 1 .4|) . By Duhamel's principle this shows 

(3.2) 6(t)= f E h {t~s)A h Q h u h , t (s)ds. 

Jo 

Using the fact that Eh(t)Ah = D t Eh(t), and Lemmas 12.11 and 12.31 we easily get 

(3.3) \\E h (t)A h Q hX \\ ^Ct-^HVQhxIl < ChH- x ' 2 \\ V X ||, for X G S h , 
and hence ^ 

\\s(t)\\ <ch 2 [ (t-sy^wvuh^Wds. 

Jo 

Here, since A^R^ — P^A, we obtain, by first applying Lemma |2.1[ 

||Vu M (s)|| < C«- 1 / 2 ||« fclt (0)|| = Cs-^WAhRnvW < Cs-^ 2 \\Av\\ = C s -^ 2 \v\ 2 , 
and hence 

\\5{t)\\ < Ch 2 f [t-s^^s^^dslv^^Ch 2 ^, 
Jo 

which completes the proof. □ 

We now consider mildly nonsmooth initial data, v £ H 1 . Here we shall need 
to assume the stability of Ph in H 1 , or ||VP/i«;|| < C|u>|i, which does not hold 
for arbitrary families of triangulations. However, a sufficient condition for such 
stability of Ph is the global quasi-uniformity of {Th}- Indeed, this assumption 
implies the inverse inequality ||Vx|| < Ch^ 1 \\x\\t which combined with the error 
bound \\RhW — w\\ < Ch\w\i, shows the desired stability of Ph- 

Theorem 3.2. Let u and Uh be the solutions of (|1.1[) and (|2 .4(1 , Then for t > 

\\u h (t) - u(t)\\ < Chh-^lvU, ifv h = P h v and \\VP h v\\ < 

Proof. Since by (|1.7p . the corresponding error estimate holds for the solution Uh of 
the standard Galerkin method (without the condition on VP^), it suffices as above 
to bound 5 = Uh — Uh- We use (|3.2|) to write 

(3.4) S(t) = {J^ + J^E h {t-s)A h Q h u Kt {s)ds = 8 1 {t)+8 2 {t). 
Using again ()3.3[) . we have, since ||Vufc it (s)|| < Cs _1 1| VP;,v|| < Cs _1 that 

\\5 2 {t)\\ < Ch 2 f (t - s)- 1 ' 2 ||Vn M (s)|| ds < Chh- 1 ' 2 Ml 

Jt/2 

Integrating by parts, we obtain 

r~ „ -lt/2 ft/ 2 

(3.5) *i(t) = E h {t - s)A h Q h u h {s) D s E h {t-s)A h Q h u h (s)ds. 
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Employing (|3.3|) . Lemmas 12.11 and 12 . 31 we now find, similarly to the above, 
(t) || < Ch 2 t-^ 2 (\\Vu h (t/2)\\ + \\VP h v\\) 

ft/2 

+ Ch 2 / (t - s)- 3 / 2 \\Vu h (s)\\ ds < C/i 2 <- 1/2 |w|i. 
Jo 

Together these estimates complete the proof. □ 

The analogous result and its proof also hold for the lumped mass method, which 
should replace the case q = 1 in [5j Theorem 3.1], since (|1.7[) does not hold for 
v h = Rhv. 

Next, we turn to the nonsmooth initial data error estimate. 

Theorem 3.3. Let u and Uh be the solutions of (|l.lj) and (|2.4|l . If (|1.17[) holds 
and Vh = PhV, then 

\\u h {t) - u{t)\\ < Cft 2 f _1 ||u||, for t > 0. 

Proof. This follows easily from the fact that for Qh satisfying (|1.17j) we have, 

(3.6) \\E h (t)A h Q h P h v\\ < Ct^WQhPhvW < ChH^Wvl for t > 0. 

This inequality is the necessary and sufficient condition for desired bound to hold 
by the following lemma, which is proved in the same way as [SJ Theorem 4.1]. □ 

Lemma 3.1. Let u and Uh be the solutions of (jl.ll) and (|2.4[) . Then 

\\u h (t) - u(t) + E h (t)A h Q h v h \\ < Ch 2 t- X \\v\\, ifv h =P h v, fort>0 7 . 

Condition (j!.17[) will be discussed in more detail in Section 0] below. Note that, 
by Lemma 12.31 without additional assumptions on the mesh, we have 

\\QhX\\<C\\VQ hX \\<Ch\\x\\, Vx&S h , 

and that the lower order error estimate of the following theorem always holds. The 
proof is the same as that of [5] Theorem 4.3]. We shall show in Section [5] that a 
0(h) bound is the best possible for general triangulation families {Th}- 

Theorem 3.4. Let u and Uh be the solutions of (|1.1|) and (|2.4[) . Then 

\\u h {t)-u(t)\\<Cht-^ 2 \\v\\, ifv h = P h v, fort>0. 

We end this section by stating optimal order estimates for the gradient of the 
error. Note that no additional assumption on {Th} is required. 

Theorem 3.5. Let u and Uh be the solutions of (|1.1[) and (|2.4[) . Then, for t > 0, 

(Ch\v\ 2 , if \\V(v h -v)\\ <Ch\v\ 2 , 

\\V(u h tt) ~ »(*))ll < \ Cht- 1 ' 2 ^ if |K - v|| < Ch\v\ u 
[Cht- 1 \\vl if v h = P h v. 

Proof. For the first two estimates it suffices, by the stability and smoothness esti- 
mates of Lemma \2. 11 to consider Vh = RhV. For this choice of the initial data the 
proofs are identical to those in [51 Theorem 3.1]. In the nonsmooth data case, the 
proof is the same as that of [5j Theorem 4.4]. □ 
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Figure 2. Left: A triangle r. Right: A patch IIf around a vertex ( 

4. Symmetric and almost symmetric triangulations 

In this section we first show that for families of triangulations {Th} that are sym- 
metric, in a sense to be defined below, assumption (|1.17|) is satisfied and therefore, 
by Theorem 13.31 the optimal order nonsmooth data error estimate holds. We shall 
then relax the symmetry requirements and consider almost symmetric families of 
triangulations, consisting of 0(h 2 ) perturbations of symmetric triangulations. In 
this case we show that (j!.17[) is satisfied with an additional logarithmic factor and, 
as a consequence, an almost optimal order nonsmooth data error estimate holds. 
Finally for the less restrictive class of piecewise almost symmetric families {Th} we 
derive a 0(h 3 ^ 2 ) order nonsmooth data error estimate. 

In addition to the quadrature error operator Qh defined in (|1.16[) we shall work 
with the symmetric operator Mh : Sh — > Sh, defined by 

(4-1) e h (iP, X ) = bP,M hX ], W>,x€S & , 

where we use the inner product 

(4-2) [1>, X ]= Yl V^M*), W^xeSV 

To determine the form of this operator, we introduce some notation. For z G Z? 
an interior vertex of Th, we define the patch IF = {Ur : r 6 Th, z £ dr}, where 
for simplicity we have assumed that r = f. Further, for z a vertex of r £ Th, we 
denote by z+ and z r _ the other two vertices of r. We then define, 

(4-3) M"* X :=~ E M(X(4)-2X(*) + X(*I)), 

for which the following holds. 

Lemma 4.1. For the operator Mh defined by (|4.1|) we have, for z G Z®, 

(4.4) Mhx(z) = M^*x with M^ X 9™en by (1431) . 
Proof. In view of (I1.16[) , we may write 

(4.5) Sh(ip,x) = (V>, JhX) ~ (*I>,X) = XI i^ J hX-^X)dx. 

reT h Jt 

For t £E Th we denote its vertices by ul,u^, and set v\ = v\, — v^, see Figure 
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[21 Writing Wj = w(yj) for a function w on t, we obtain, after simple calculations, 
(4-6) [ i>JhXdx = j^Yl ^ 22 ^ + 7 Xj-i + 7 Xi+i), 



and 



Thus 



/| T | 3 
^Xdx = — !pj( 2 Xj + Xj-i + Xj+i) 

11 3=1 



3 



{^JhX-^x) dx = —^j2^,ipj(Xj+i - ZXj+Xj-i)- 



54 . 

3=1 

Summation over r £%,, (|4.ip and (|4.5[) show 

[V,M^ X ] = J! <P(z) M h*X, W>,x G S^- 

This implies (|4.4I) and thus completes the proof. □ 

We say that Th is symmetric at z £ Z®, if the corresponding patch Tl z is sym- 
metric around z, in the sense that if x £ II 2 , then z — (x — z) = 2z — x G II 2 . We 
say that Th is symmetric if it is symmetric at each zG^JJ. The patch IIf in Figure 
[5] is nonsymmetric with respect to £o; whereas triangulations which are built up of 
either of the patches shown in Figure [3] are symmetric. Symmetric triangulations 
exist only for special domains, such as parallelograms, but not for general polygonal 
domains. 

We now show the sufficiency of symmetry of {Th} for condition (|1.17[) for the 
operator Qh, and hence, by Theorem 13. 31 for the nonsmooth data error estimate. 

Theorem 4.1. If the family {Th} is symmetric, then (| 1 . 1 T[) holds. 

Proof. The proof, by duality, follows that of [5j Theorem 5.1]. For given \ G Sh we 
define ip — <p x £ H 1 as the solution of the Dirichlet problem — Ap = % in fi, <p = 
on dil. Since ft is convex, we have ip £ H 2 and \<p\2 < C||xll- With Ih the finite 
element interpolation operator into Sh, we have, for any tp £ Sh, 

(4.7) \\Qhm\ = sup — — — = sup — 

x eS h llxll x es h llxll 

\(VQ h ^,V{ V -I h p))\ KVQ^W^_ 

-- sup 1 1 1 1 ' II II — \ J- J- • 

xes h llxll x&s h llxll 
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By the obvious error estimate for Ih and Lemma l2.3[ with p = 0, we find 

(4.8) \I\<Ch sup l|V y' IM2 <Cfe 2 |H|. 

xes h llxll 

To estimate II, we employ (|1.15[) and (14. 1|) to rewrite the numerator in the form 

(4.9) {VQ h ip,VI h <p) = e h {^,I hV >) = [i/j,M h I h (p]. 

To bound M^Ihip, we consider an arbitrary vertex z = (q G Z®. Let H^ be 
the corresponding patch of 7/j, with vertices {CjljLi? numbered counter-clockwise, 
with Q+K = Q for all j. Also denote by {rj}f =ll the triangles of 7^ in U^ , with 
Tj having vertices Co> Cj, O+i? an d set r o = T if (see Figure [2]). Then Lemma [4. II 
implies 

1 K 

(4.10) M h lMCo) = M^°I W = --^^(^(0) - vj(Co)), 

3=1 

with Wj = |tj_i| + |rj|. By assumption, the patch Il^ is symmetric and hence, 
by (I4.10p . we can express MhIh<p((o) as a linear combination of terms of the form 
(p(Cj) — 2ip(Co) + <p{Cj)-> where Co is the midpoint of the vertices Q and Cj of U^ . 
Hence M^Ih^iCo) = for ip linear in IT^ and, as in [5], we may apply the Bramble- 
Hilbert lemma to obtain 

(4.11) \M h lM(a)\ < Ch 2 \n Co \ 1/2 MH H n (a ) < C7i 3 |M| ff3( n Co ). 
Employing this estimate for all patches Ll z of Th, we obtain, for any i/ifS/,, 

(4 12) \[1>,M h I h <p]\ < Ch 3 Mz)\ \\<p\\n>(p M) < Ch 2 \\iP\\ \<p\ 2 < Ch 2 U\\ || X ||. 

Hence, in view of (|4~7| and (jUSJ, we obtain |J7| < C/z. 2 1 1 V 7 1 1 - Together with (|4~8]) 
this completes the proof. □ 

We now want to slightly weaken the assumption about symmetry. We say that 
a family of triangulations {Th} is almost symmetric if each Th is a perturbation by 
0(h 2 ) of a symmetric triangulation, uniformly in h, in the sense that with each 
patch n z of Th there is an associated symmetric patch from which n z is obtained 
by moving each of its vertices by 0(h 2 ). Such triangulations exist for any convex 
quadrilateral, cf. Figure 01 We note that various special triangulations have been 
used in the past for obtaining higher order accuracy for the gradient of the finite 
element solution (super-convergent rates of 0{h 2 ) or 0(h 2 £h j), see, e.g., [71 IT!] [16] . 
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For example, the strongly regular triangulations from requiring that any two 
adjacent triangles form almost a parallelogram (a deviation of a parallelogram by 
0{h 2 )), are almost symmetric meshes in our terminology. We shall show that, in 
this case, we have almost optimal order convergence for nonsmooth initial data. 

Theorem 4.2. If the family {Tk} is almost symmetric, then 

(4.13) HQ^II < Ch 2 t][ 2 HVH, V^e5 k , where £ h = l+\logh\. 
Hence, for the solution of (|1.13|) . with Vh = PhV, we have 

(4.14) \\u h (t) -u{t)\\ < Ch 2 i]l 2 t- X \\v\l fort>0. 

In the proof we shall need the following Sobolev type inequality, where the | ■ \ H k 
denote seminorms with only the derivatives of highest order k. 

Lemma 4.2. Let B be a fixed bounded domain, satisfying the cone property. Then 
we have, for < e < 1, 

\tp(z') - tp(z)\ _ 1/2 , . 2 

z,z'eB, z'=£z \ z — z \ 

Proof. We find from [TJ pp. 109-110], for e small, with C independent of e, 

(4.15) sup ^-.ff 1 <C||Vy.|| MB) , withp = 2/e, V<peW}(B). 

z,z'£B, z'jtz \ z ~ z \ 

We shall also apply the Sobolev inequality, with explicit dependence on p, 

(4.16) Ml p ( B ) < Cy/ 2 |MU 1( B), forp< TO , V^EH l {B). 

For ip G Hq(B) a proof was sketched in |T5} Lemma 6.4]. For the general case of 
(p G H 1 (B), we make a bounded extension of ip from H 1 (B) to Hq(B), with B C B, 
cf., [3 IV] and apply (|4"?TB1) to H l (B) to complete the proof. 
Employing (|4.16l) yields 

||V^|| ip(B) < C^ 2 {\ip\ B x {B) + W\hhb)), g H 2 (B). 

Combining this with (|4.15p . using p 1 / 2 = (2/e) 1 / 2 , completes the proof. □ 



Proof of Theorem \4-2\ The proof proceeds as that of Theorem 14.11 starting with 
(|4.7p and noting that the bound (|4.8j) for / remains valid. In order to bound //, 
we follow the steps above, but now, instead of (|4.11[) . we show 

(4-17) \M h I hl p(Co)\ < Ch% /2 M H i (U(o y 

Using (|4~T7)) as (|4~TT|) in (|4T2)) . we find 

(4.18) \[4>,M h I h <p]\ <Ch 2 i)! 2 U\\ llxll, V^xeSi, 

and hence \II\ < Ch 2 £l/ 2 \\il>\\- Together with (|4.8|) . this completes the proof of 
(|4.13p . The error estimate (|4. 14|) now follows from Lemma T3. II and 

\\E h (t)A h Q h P h v\\ < Ct- 1 \\Q h Phv\\ < Ch 2 l]l 2 t- l \\v\\, for * > 0. 

It remains to show (|4.17[) . Let LT^ be the symmetric patch associated with Ii^ a 
by the definition of almost symmetric. After a preliminary translation of 11^/ by 
0(h 2 ), we may assume that £q = £o- Further, without loss of generality, we may 
assume that U^ C H^ . In fact, if this is not the case originally, it will be satisfied 
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by shrinking IL^ by a suitable factor 1 — ch 2 with c > 0. Starting with IL^ we 
may now move the vertices one by one by 0(h 2 ) to obtain IL^ in a finite number 
of steps, through a sequence of intermediate patches H^ C n^ . 
Applying (|4.10[) we will show that for each of these 

(4.19) \M^ <0 I h <p\ < C e h 3 - e \\<p\\ H 2 {nCo) , where C e = Ce^ 2 , e > 0, 

which implies (|4.17l) . by taking e = l^ 1 and YVq = II^ . 

Since (14. 19|) holds for the symmetric patch 11^ , by (|4.11[) . it remains to show 
that if it holds for a given patch H^ then it also holds for the next patch in the 
sequence. Assuming thus that (|4.19l) holds for H^ , we consider the effect of moving 
one of its vertices, £2, say, to £2) with \( 2 — ( 2 \ = 0{h 2 ). 

Applying Lemma l4~2l to the function ip(h-), with B suitable, we obtain 



(4.20) sup 1^1 ff 1 < C e h-^{\ip\m { n i0 ) + %|^(n Co) ) 



< c e h- 1+c \ 



Moving only the vertex £2 in n< changes only the triangles t± and t 2 and thus the 
terms corresponding to j = 1, 2,3 in (|4.10[) . 

Letting t[ and t 2 be the new triangles, the change in the term with j = 1 is then 
bounded, since ||r(| — |ti|| < Ch 3 , by 

WCO-^Co)! 



K-«i)(v>(Ci) - v(Co))| < C\\t[\ - Inl^ 1 - 

< C £ /l 3 ||^|| ff 2 (II( . o ) 



Id - Col 1 " 



and thus by the right hand side of (|4.19p . The change in the term with j = 3 is 
bounded in the same way. For j — 2 the change is bounded by the modulus of 

W2(V(^)-V(0)))-W2(V«2)-V(CO)) 

= K - wa)(p(Ca) - p(Co)) + <4(<P(Q ~ 
The first term on the right is bounded as the terms with j = 1,3, and the second 
is bounded, using (|4.20p . since C2 — C2 1 < Ch 2 , in the following way 

\<SMQ - v(Ca))| < C £ h 2 \( 2 - C2\ 1 - e h- 1+e \\i P \\ H 2 (Il<o) < C e h 3 - e M HHn<o) . 

This shows that f|4. 19[) remains valid after moving £2 , which concludes the proof. □ 

More generally, we shall consider families of piecewise almost symmetric trian- 
gulations {Th}, in which Q is partitioned into a fixed set of subdomains {Qf-}k=i, 
and each of these is supplied with an almost symmetric family {7ft, (fife)} so that 
Th = U^ =1 7ft,(r2fe). Such families may be constructed for any convex polygonal 
domain, cf. Figure [4j by successively refining an initial coarse mesh, a procedure 
routinely used in computational practice. For such meshes we show the following 
result. 

Theorem 4.3. If the family {Th} is piecewise almost symmetric, then 

(4.21) \\Q h iP\\<Ch 3 / 2 U\\, WeS h . 
Hence, for the solution of (|2.4[) with — PhV, we have 

(4.22) \\u h {t) -u{t)\\ < C7i 3/2 i -1 |M|, fort>0. 
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Proof. Following again the steps in the proof of Theorem 14. 1[ we note that (|4.8j) 
still holds, and it remains to bound II. For each internal vertex (o of one of the 
Th{Qk), the corresponding patch YIq is a 0(h 2 ) perturbation of a symmetric patch, 
and thus (I4.17P holds. For Co S a vertex on the boundary of two of the Th(&k) 
we see that by (|4.10|) 

\M hX (Co)\ < Ch 3 max |Vx(ar)| < Ch 2 \\ V X \\l 2 (u <0 ), 

X ^ Co 

and by the use of approximation properties of the interpolation operator Ih we get 

(4.23) \M h lMCo)\ < Ch 2 \\I w \\ m{n(o) < Ch 2 (\\<p\\ m(I1(o) + %|^(n <0 )). 
Using ([1117} and P~2"3")) as earlier (|4"lTt in P~12")) . we conclude 

\ty : M h i h ^\<ch 2 il /2 u\\\p\ 2 + chU\\y\\ HHns) , 

where fig is a strip of width 0(h) around the interface between the subdomains £lk 
of fi. Using now the inequality |M|jji(n s ) < Ch 1 / 2 \\ip\\ H *{n) < C7i 1/2 ||xll, we get 

(4.24) \bP,M h I h <p]\ < Ch 3 / 2 \\yj\\ Hxll, W>,x G 5 h5 

and hence |77| < C7i 3/2 1|^||. Together with (|4.8[) . this completes the proof of (|4.21l) . 
The error estimate (|4.22[) now follows by Lemma 13.11 and 

||^ fc («)A h Q h i\«|| < Ct^WQhPuvW < Ch*tH- 1 \\v\\, for <> 0. □ 

We remark that the operator Mh used here, modulo a constant factor, is the 
same as the operator in [5] . The arguments in the proofs of Theorems 14.21 and 
14.31 therefore show that the following result holds for the lumped mass method. 

Corollary 4.1. Assume that {Th} is almost or piecewise almost symmetric. Then 
the nonsmooth data error estimates for the lumped mass method, corresponding to 
(|4.13p and (|4.21[) . respectively, hold. 

We finish this section by remarking that, in one space dimension, the full 0(h 2 ) 
Li norm bound (|1.17p for Qh holds also for almost symmetric partitions, without a 
logarithmic factor. Let fl = (0, 1) be partitioned by = xq < x± < ■ ■ ■ < x^ h +i 
1 . Denote now Th = {tj}^-^ , with Tj = [:Ei_i,£j], and let Sh be the set of the 
continuous piecewise linear functions over Th, vanishing at x = 0, 1. We set hi = 
Xi—Xi^i and h = max^ hi. The control volumes are Vj = (xi~ hi/2, Xi~\-hi+i/2) and 
Jhip{x) = ip(xi) for x G Vi- We say that Th is almost symmetric if — hi\ < Ch 2 
for all i. 

Simple calculations show, with (XiVO — J xV' dx and (x,^) = (x,Jhip), for 

x,V> e S h , 

i=i 

where Wi = w(xi) for a function w on 17, and the one-dimensional version of (|4.10[) 
at Xi becomes 

Mhlh<p(xi) = --^(hi + i(ipi+i - tpi) + hiQpi-i - tpi)), i = l,..., N h . 
The crucial step to prove (|1.17j) is then to show an analogue of (|4.1ip . in this case 

(4.25) \M h I h <p{ Xi )\ < Ch^ 2 y\\ H 2 {n ^ }l i = l,...,N h , with H,., =r i Ur m , 
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Figure 5. Left: A nonsymmetric mesh. Right: A nonsymmetric 
patch n^ , around £q. 



from which (|1.17| follows as earlier. Using the Taylor formula 

<p(x) = (p(xi) + (x - Xi)ip'{xi) + / (x - y)(p"(y)dy, 

J Xi 

we find easily 

M h I h <p(x 4 ) - ~(h 2 i+1 - h^cp'ixi) + 0(h 5 / 2 H^'IU,^)), i = l,...,N h . 
By the almost symmetry, \hj +1 — hf \ < C/i 3 and by the Sobolev type inequality 

\<p'(xi)\ < ^ 1/2 (||^|| L2( n Xi) + hy\\ L2i u Xi) ) < Ch-^Mn^n^, 
for i = 1, . . . , JV/,, we now conclude that (|4.25l) holds. 

5. Examples of nonoptimal nonsmooth initial data estimates 

In this section we present two examples where the necessary and sufficiency 
condition (J3T6J) for an optimal 0(h 2 ) nonsmooth data error estimate for t > is 
not satisfied. In the first example we construct a family of nonsymmetric meshes 
{Th} for which the norm on the left hand side of (|3.6p is bounded below by ch, 
thus showing that the first order error bound of Theorem 13.51 is the best possible. 
In the second example we exhibit a piecewise symmetric mesh for which this norm 
is bounded below by c/i 3 ^ 2 , implying that the error estimate of Theorem 14.31 is best 
possible. 

In our first example we choose Q — (0, 1) x (0, 1) and introduce a quasiuniform 
family of triangulations {Th} of ft as follows. Let N be a positive integer divisible 
by 4, h = 4/(3iV), x = 0, and set, for j = 1, . . . , N and m = 0, 1, . . . , M = |iV, 



(5.1) 



Xj-l 




for j odd, 
for j even, 



and 



mh. 



We split the rectangle (xj,Xj+i) x (y m ,ym+i) into two triangles by connecting the 
nodes (xj, y m ) and (xj+i, y m -i), see Figure [5l This defines a triangulation Th that 
is not symmetric at any vertex. 

Let now £o = (^2j, ym), Co S Z^, and let II^ be the corresponding nonsymmetric 
patch shown in Figure with vertices {Cj}?=i- Let tj be the triangle in Hq with 
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vertices Co, (j, Cj+i, where (V = Ci- We then have |Tj| = \h 2 , for j = 1,2,3, 
and \Tj\ = 5ft 2 , f° r j — 4,5,6. Thus, using (|4.10[) . for %j) £ Sh, we obtain with 

= i>(Q)i 

1 ^ 1 

M fc ^(Co) = - ~ £>iW»i - iM = -T7-r(W>i +i>4- 2^o) 
(5.2) 54 j=i 54 4 V 

+ 2(^2 - V'o) + 2(^3 - V>o) + 4(^5 - V>o) + 4(^6 - V>o)) ■ 
Because V-0 is piecewise constant over II^ , we easily see that (|5.2p implies 
(5-3) \M,MCo)\ <C^ 2 ||V^|| i2(n<0 ), W>eS h . 

For a smooth function <p we have, by Taylor expansion, 

- ^(Co) = V^(Co) ■ (0 ~ Co) + O(^), 

where £j is considered as a vector with components its Cartesian coordinates and 
the dot denotes the Euclidean inner product in W 2 . Employing this in (|5.2[) . we 
find, after a simple calculation, 

(5.4) M h I hV >(Co) - ^V^(Co) • (3, -1) + 0(ft 4 ). 

Let 4>%(x, y) = 2 sin(7ra;) sin(-7ry) be the eigenfunction of — A, corresponding to the 
smallest eigenvalue Ai = 2n 2 . We then easily find that V</>i(l/4, 1/4) • (3, —1) = 2tt. 
Hence, there exists a square V = [1/4 — d, 1/4 + d] 2 , with < d < 1/4, such that 

(5.5) V<h.{z) ■ (3,-1) > 1, VzeV. 

Letting now for z £ Z\ PiV we then have that Mhlh4>i{ z ) > eft 3 , c > 0, for ft, small. 
We shall prove the following proposition. 

Proposition 5.1. Let Th be defined by (|5.ip . Vh = {z = (x2j,y m ) £ V} and 
consider the initial value problem (|2.4[) with Vh — ^2 z ^-p h where $ 2 £ Sh is the 
nodal basis function of Sh at z. Then we have, for h small, 

\\E h (t)A h Q h v h \\ > c(t)h\\v h \\, with c{t) > 0, for t > 0. 

Proof. Letting Xj and be the eigenvalues and eigenfunctions of — A^, and using 
Parseval's relation in Sh, equipped with (•,•), we have 

(5.6) \\\E h (t)A h QhV h \\\ 2 =J2 e ~ 2tl H A hQhV h ,$} 2 > e- 2tx *(A h QhVh,rt) 2 . 

3 = 1 

Combining (f273|) . (fl~15|) and (|4~T|) . we find 

(5.7) - (A h Q h v h ,ip) = (VQ h v h ,Vil)) = e h (v h ,ip) = [vh,M h ip], Vi/> G S h . 

Note now that for z £ Vh, the corresponding patch H z has the same form as the 
patch il£ considered above. Thus employing (|5.3I) for Co = z we get, for ip £ Sh, 

(5.8) \[v h ,M h n < ^ J2 IWI < Cft||V^||, 

zev h zev h 

where in the last inequality we have used the fact that the number of points in Vh 
is 0{N 2 ) = 0{h~ 2 ). We recall from [W\ that 

\\4>i - <Pi\\m = 0(h) and Ai -> Ai, as ft ^ 0, 




Figure 6. A piecewise symmetric mesh. 

and, since obviously ||0i — Ii^iWh 1 = 0(h), (|5.8p with tp = <ffi — Ih(f>i, gives 

(5.9) \[v hl M h {j ) h l -I h d> 1 )]\<Ch\\V{j> h l -I h (t )1 )\\ <Ch 2 . 

For every z G V%, (|5.5I) holds, and thus, using (I5.4I) with = and Co = z , we 
obtain, for h small, since the number of vertices in Vh is bounded below by cN 2 , 

[v h ,M h Ih<h.] = M hhM z ) > ch3]S[2 = ch > with c > °- 
Combining this with (|5.9p . we obtain, for h small, 

[v h , M h $] > [v h , M h I h 4>x] - \[v h ,M h {4>\ - I h <t>i)]\ >ch- Ch 2 > ch, with c > 0. 
Since \\\v h \\\ = O(l), (f5T6|) and (|5J| now show 

\\\E h (t)A h Q h v h \\\ >e-^[v h ,M h ^] > c(*)/i |||u fc |||, for t > 0. 
Since ||| • ||| and || ■ || are equivalent norms, the proof is complete. □ 

It follows from Proposition 15 . II and Lemma |3. II that the highest order of conver- 
gence that can hold, uniformly for all u £ L2, and for any family of triangulations 
{Th}, is O(h), i.e., Theorem l3.4l is best possible, in this case. 

We now turn to our second example, in which {Th} is a piecewise symmetric 
family. Let again fl = (0, 1) x (0, 1) and consider a triangulation Th of fl, where the 
nodes (xj, y m ) are given as follows. With J a positive integer, let N — 7 J, M = 4 J 
and h = 1/(4 J), and set for j = 0, . . . , N and m = 0, . . . , M, 

(5.10) Xj = [ jh ; for0<i<J, 

V ' 3 \l/4+(i- J)h/2, for J <j<N, 

see Figure |U This time we consider the set of vertices in V with x = 1/4 and prove 
the following proposition. 

Proposition 5.2. Let Th be defined by (|5.10[) and T 5 ^ = {z = (xj,y m ) G "P}. For 
the initial value problem (|2.4[) , = X^^ep' ®z> where $ z G ^ is the nodal 

basis function of Sh at z, we have, for h small, 

\\E h (t)A h Q h v h \\ > c(t)h 3 / 2 \\v h \\, with c(t) > 0, for t > 0. 
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Proof. Again, using (|5.6j) and (|5.7[) . we have, 

(5.11) \\\E h (t)A h Q h v h \\\ 2 > e- 2t ^[v h ,M h 4>1} 2 . 

For z G T" h , the corresponding patch II z has the same form as the patch n^ 
considered above, see Figure [5] (right). Thus employing (|5.3I) for (o — z and taking 
into account that the number of vertices in V' h is 0(N) we now obtain, for i\> 6 Sh, 

\[vh,M h i/>]\ < \[$*,M h il)]\ = \ M hi>(z)\ < Ch 3 / 2 \\Vip\\. 
zev' h zer' h 

Similarly to (15.91) this now shows 

(5.12) \[v h ,M h $l-I h <t>x)]\ <Ch 5 / 2 , 
and, again using (I5.4j) , for h small, 

MfcJ h 0i] = M hh<f>i{z) > ch 3 J = ch 2 , with c> 0. 

Combined with f|5 . 1 2[) this gives, for h small, 

(5.13) [v h , M h $] > ch 2 - Ch 5/2 > ch 2 , with c> 0. 
Since H^H = 0(h}/ 2 ) we obtain from (|5TTTl and (|5~13)) 

\\\E h (t)A h Q h v h \\\>c(t)h 2 >c(t)h 3/2 \\\v h \\\, fori>0. □ 

It follows from Proposition 15.21 and Lemma 13.11 that the highest order of con- 
vergence that can hold, uniformly for all v £ L2, and for all piecewise symmetric 
families {%,}, is 0(/i 3/2 ), i.e., Theorem S3] is best possible in this regard. 

Remark 5.1. Since Mh is proportional to the operator used in [5], the ar- 
guments in this section also apply to the lumped mass method. In particular, the 
analogue of Proposition I5.il then shows that the first order nonsmooth data esti- 
mate for t > of [SJ Theorem 4.3] is best possible for general triangulations {Th\. 
Further, the 0(h 3 ^ 2 ) estimate stated in Corollary \4-l\ is best possible for piecewise 
almost symmetric triangulations. Our examples here may be thought of as gener- 
alizations to two space dimensions of the one-dimensional counter-examples in [5J 
Section 7] . 

6. Some fully discrete schemes 

In this section we discuss briefly the generalization of our above results for the 
spatially semidiscrete finite volume method to some basic fully discrete schemes, 
namely the backward Euler and Crank-Nicolson methods. 

With k > 0, t n — nk, n = 0, 1, . . . , the backward Euler finite volume method 
approximates u(t n ) by U n £ S h for n > such that, with dU n = (U n - U n - 1 )/k, 

(8U n ,x) + (VU n ,V X )=0, VxeS ft ., forn>l, with U° = v h , 

or, 

(6.1) dU n - A h U n = 0, forn>l, with U° = v h . 
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Introducing the discrete solution operator Eun = (J — kAh) 1 we may write 
U n = EkhU"^ 1 = E'£ h U°, n > 1. Using eigenfunction expansion and Parseval's 
relation, we obtain, analogously to JTSJ Chapter 7] , the stability property 

(6.2) \\V p E2hX\\ < C\\V P X\\, Vxe5 h) for p = 0,1. 

The estimates that follow and their proofs are analogous to those for the lumped 
mass method derived in [5], since the operators Eh(t), Ah and Qh, defined in Section 
[2j have properties analogous to those of the corresponding operators for the lumped 
mass method. For simplicity we will only sketch the proof of Theorem 16.11 

We shall use the following abstract lemma shown in [5], in the case H — Sh, 
normed by ||| • |||. and with A = — Ah- 

Lemma 6.1. Let A be a linear, self adjoint, positive definite operator in a Hilbert 
space T~i, with compact inverse, let u = u(t) be the solution of 

u' + Au = 0, for t > 0, with u(0) = v, 

and let U = {U™}^L be defined by 

<9U" + AU"=0, forn>l, withV a =w. 

Then, for p = 0, 1, — 1 < q < 3, with p + q > 0, we have 

||AP/ 2 (U"-u(i„))|| < CTrf-^-^IIA^+^vll, forn> 1. 

The error estimates of the following theorem for (|6.1I) are of optimal order under 
the same assumptions as in Section [3] 

Theorem 6.1. Let u and U be the solutions of (jl.ip and (|6.ip . Then, for n > 1, 

(C(h 2 + fc)Ma, if \\v h -v\\ <Ch 2 \v\ 2 , 

\\U n -u(t n )\\ < lc(h 2 +k)tn 1/2 \v\ 1 , if Vh = P h v and \\\7P h v\\ < C\v\ u 
[C(h 2 +k)t- 1 \\v\\, if Vh = P h v and $UJi) holds. 

Proof. Analogously to the proof of [SJ Theorem 8.1], we split the error as 

U n - u{t n ) = (U n - u h {t n )) + (u h (t„) - u(t n )) =p n + r) n . 

By Theorems I3.1[ 13.21 and 13. 3[ r\ n is bounded as required. In order to bound 
f3 n = (E£ h — Eh{t n ))vh in the smooth data case, it suffices, using the stability 
estimates (|6.2p and Lemma \2. 11 to consider Vh — RhV. We obtain by Lemma \6. 11 
with A = A h = -A h , and q = 2, 1, 0, 

|||/3„||| = |||^ n -u h (i„)||| <Ckt-^-^\\\A q h /2 v h \\\ <Ckt-^-^\v\ q , 
where for q = 2, the last inequality follows from 

\\\A h R h v\\\ 2 = (\7R h v,VA h R h v) = (Vv,VA h R h v) = -{Av,A h R h v), 

for q = 1 from \ \\Al /2 P h v\\\ = \\VP h v\\ < C\v\i and for q = from |||J\»||| < 
C\\v\\. " □ 

Also for the lumped mass method the analogous result in the mildly nonsmooth 
data case v £ H 1 holds, and should replace the result for q = 1 in [5J Theorem 8.1], 
cf. the remark after Theorem 13.21 

Recall that Qh satisfies (|1.17j) if {Th} is symmetric. For almost symmetric or 
piecewise almost symmetric {7h} we obtain correspondingly the following non- 
smooth initial data error estimates employing (|4.14|) and (|4.22|1 . 
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Theorem 6.2. Let u and U be the solutions of (11. 1|) and (|6.ip . with Vh = PhV. 
Then, for n > 1, 



\\U n -u(t n )\\ < 



C{h 2 i h ^ 2 + /c)t n 1 ||w||, if {Th} is almost symmetric, 
C(h 3 / 2 + fe)*n 1 |l u llj if {Th} is piecewise almost symmetric. 

For the gradient of the error we may prove as in [5j Theorem 8.2], the following 
smooth and nonsmooth data error estimates, without additional assumptions on 
Th- For smooth initial data we assumed in [5] that Vh — Rh v , but the more general 
choices of Vh are permitted by the stability estimates (|6.2[) and Lemma 12.11 

Theorem 6.3. Let u and U be the solutions of fll.lj) and (|6.ip . Then, for n > 1, 

III (nm -\c(ht^ + kK 3/2 )\\v\\, if v h = P h v. 

We now turn to the Crank-Nicolson method, defined by 
(6.3) dU n -A h U n -\ =0, forra>l, with [7° = %, C/"-^ = ±(J7 n + t/"" 1 ). 

Denoting again the discrete solution operator by E kh = (I + |fcA/ l )(7 — ifcA/j) -1 

we may write ?7 n = EkhU" 1 ^ 1 = EJ! h U , n > 1. Using eigenfunction expansion and 
Parseval's relation, we find that (16.21) also holds for this method. 



The Crank-Nicolson method does not have as advantageous smoothing proper- 
ties as the backward Euler method, which is reflected in the fact that the following 
analogue of Lemma 16.14 shown in [5j Lemma 8.2], does not allow q = 0. 

Lemma 6.2. Let A and u(t) be as in Lemma \6.1\ and let U n satisfy 
BV n + AU""5 =0, forn>\, withU°=v. 

Then 

|| A p/2(un _ u (t n ))|| < Ck 2 t- (2 - q) \\A p / 2+q v\\, forn>l, p = 0, 1, q = l,2. 

This time optimal order estimates for the error in L 2 an d in H 1 , hold uniformly 
down to t = 0, if v G H 4 and v g i/ 5 , respectively. The proofs are analogous to 
those of [5l Theorems 8.3 and 8.4], where we assumed Vh = RhV- Again the stability 
estimates (|6.2I) and Lemma |2~T1 permit the more general choices for Vh- 

Theorem 6.4. Let u and U be the solutions of (|1.1|) and ()6.3|) . Then, with q = 1, 2, 
we /lave, /or n > I, 

\\U n - u(t n )\\ < C(h 2 + k 2 t- {2 -^)\v\ 2q , if \\v h - v\\ < Ch 2 \v\ 2 

\\V(U n - u(t n ))\\ < C(h + k 2 t-^)\v\ 2q+1 , if \\V(v h -v)\\< Ch\v\ 2 . 

For optimal order convergence for initial data only in L 2 , one may modify the 
Crank-Nicolson scheme by taking the first two steps by the backward Euler method, 
which has a smoothing effect. We may show then the following result, analogously 
to that of Theorem 8.5], with the obvious modifications for almost symmetric 
and piecewise almost symmetric families {Th}- 

Theorem 6.5. Let u be the solution of (II. ip and U n that of ()6.1j) . for n = 1,2, 
and of (|6.3j) . for n > 3, with Vh — PhV and assume (|1.17[) holds. Then we have 

\\U n - u{t n )\\ < Cihh- 1 + k% 2 )\\v\\, for n>l. 
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7. Problems with More General Elliptic Operators 

This final section is devoted to the extension of our earlier results to the more gen- 
eral problem (I1.19[) , and we recall that we shall consider the finite volume method 
(|1.26p where the bilinear form «>,(■, •) is defined by (|1.25l) . Our error analysis is 
again based on estimates for the standard Galcrkin finite element method, in this 
case defined by (|1 .20|) and (|1.2ip . It is well known that for this method the stability 
and smoothing estimates (|2.2j) hold as do the error estimates ()1.5j) - (11.7|) . where the 
norms | • \ q are defined analogously to the norms (11.21) . using the eigenvalues and 
eigenfunctions of A. 

We introduce the discrete elliptic operator Ah : Sh — > Sh by 

(7-1) (A h ip,x) = a h (ip, J h x), Vx,Ve5ft, 

which is symmetric and positive definite with respect to the inner product (•, •) by 
(|1.24l) . since (/3ip, JhX) is symmetric, positive semidefinite on Sh- This follows from 
the fact that J X^hip dx is symmetric by (|4.6p and ft is constant and nonnegative 
in each r of 7/,. We may then rewrite (|1.26p as 

(7.2) Uh,t + A h u h = 0, for t > 0, with u h (0) = v h , 

and the solution is given by Uh(t) = Eh(t)vh, where Eh(t) — e~ Ah * is defined as in 
(|2.5p . with {Ay} and {</>^} the eigenvalues and eigenfunctions of Ah, orthonormal 
with respect to (•,•). 

Note that a slightly different finite volume element method for (|1.19p has been 
considered in |12j . This method differs in the discretization of the lower order term, 
using the bilinear a^O, •) defined by 

BfcM JhX) = (svv. v x) + WJhi>, JhX), v^, X G s h . 

For this method analogous results to Theorems 17.1 l - [?T3l hold . 

Following our error analysis in the previous sections we introduce 5 = Uh — Uh 
and split the error into Uh — u = S + {uh — u), where uu — u and W(v,h — u) are 
estimated by the analogues of (|1.5p - (|1.7p . It therefore suffices to derive estimates 
for 5, which satisfies, for t > 0, 

(7.3) (8 h ,t,x) + ah(8,JhX) = -£h(uh,t,x)-£h(uh,X), V% £ Sh, with 6(0) = 0, 
where e/ l (-, •) is given by and £/,(•,•) is defined by 

(7-4) £h(ip,X) = ah(ip,JhX)-a(' l P,x), V^xeS^. 

Now let Qh : Sh ^ Sh and Qh ■ Sh — > Sh be the quadrature error operators given 

by 

(7.5) a h (Qh^,JhX)=£h(^,x) and a h (Q h ip,J h x)=£h(^,x), V-0, X ^ 
Using (|7.ip . the equation (|7.3p for 6 can then be written in operator form as 
5 t + A h 6 = -A h Q h u htt -A h Q h u h , for t > 0, with 5(0) = 0. 

This problem is similar to (|3.ip . except that the operator —Ah is replaced by Ah 
and that on the right hand side we have an additional term resulting from the 
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approximation of the bilinear form a(-, •). By Duhamel's principle we have 
5(f) = - / E h (t - s)A h Q h u Kt {s) ds 



(7.6) 



E h {t - s)A h Q h u h {s) ds =: 5(t) + 6(t), for t > 0. 



To estimate 5 it therefore suffices to bound S and <5. For this we need some 
auxiliary results, which are discussed below. 

Lemma 7.1. Let a, /3 G C 2 . For the error functional Sh, defined by (17.41) . we have 

\e h {i>,X)\ < C^+1VVII l|V p X ||, W> )X eS fe) With p,q = 0,1. 
Proof. In view of (|7.4[) . we may write 

£/>Wsx) = ((5 - a)VV>, V%) + 0ip,JhX) ~ W,x)- 
We then split e^^x) as a sum of integrals over r G Th- Since a = a(z r ), we see 
J T (f — f(z T ))dx = for linear functions /, and hence 

(7.7) \J(f-f(z T )dx\<Ch 2 T \r\\\f\y, ^ feC 2 , 

with h T the maximal side length of r. Therefore, using this and the fact that 
Vip ■ V% is constant in t, wc get 

| J {a - a)W ■ V X dx\ < Ch 2 T \\a\y J |V^ • V X | dx < Ch 2 T \\V^\\ L2(T) \\V X h 2 (r)- 
Employing an inverse inequality locally and summing over r £ Th, we obtain 

(7.8) |((3 - a)W, Vx)| < C/i^H Wll ||V P xll- 
In a similar manner we estimate the zero order term. Obviously, 

(7.9) (W,i/ i x)-W,x) = £/ i (#,x) + P-W,x)- 

Using Lemma 12.21 we can bound the first term on the right-hand side of (|7.9j) , as 
desired. We then split the second term, in the following way 

0-f3)iP X dx= [ - P)ty x ){z T )dx 



(7.10) 

+ J (fi-0)(il>x-W>x)(zr))dx=:I + II. 

Employing (|7.7p we easily get 

\I\ < Ch 2 T \\f3\y\r\\^x)(zr)\ ^ ChUr]- 1 ] f ^dx\ I J X dx\ 

<Ch 2 U\\ L2{T) \\x\\L 2{ r), 

and since \j3 — /3| < Ch T \\/3\\ c i in t, 

|77| < Ch 2 T J (\Vtp X \ + \^^ X \)dx 

< Ch 2 (\\Vyj\\ L2{T) \\x\\ L2 (r) + ||^lk(r)l|VxlU a (r))- 

Combining the bounds for / and // with (|7.10j) . using an inverse inequality locally, 
summing over r e Th and using (|7.8j) . we conclude the proof. □ 
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For the solution operator Eh(t) = e Ah 4 of (|7.2|) . one shows, as in Lemma \2.1\ 
the following smoothing property. 

Lemma 7.2. For Eh, the solution operator of (|7.2|) . we have, for Vh € Sh and 

t > 0, 

\\V p D e t E h (t)v h \\<Ct- e -^/ 2 \\V^v h ie>0, p,? = 0,l, 2^ + p>g. 

Further, following the steps in the proof of Lemma 12.31 we can get easily the 
following estimate 

Lemma 7.3. Let Ah, Qh andQh be the operators defined by (|7.ip and l|7.5[) . Then 

\\WQhx\\ + h\\AhQhx\\<ChP +1 \\VP x \\; VxeSh, for p = 0,1, 
and the same bounds hold if we replace Qh by Qh. 

Proof. Using the fact that a,h(x, JhX) > c ll^xl| 2 ; f° r X G $h, <|7.5|) and Lemma [2~2l 
with ip = QhXi we obtain for p = 0, 1, 

c||VQ,xl| 2 < o A (Q AX) JfcQfcX) = O^X) < C7^ +1 ||V p xll \\VQhXl 

which bounds QhX as desired. By the definition of Ah and Lemma 12.21 with q = 0, 
we also get for p = 0, 1, 

III^QhxIII 2 = ehfaAhQhx) < ChP\\^p x \\ \\A h Qhx\\. 

Since the norms ||| • ||| and || • || are equivalent on Sh, this shows the bound stated. 

To prove the corresponding bounds for Qh, analogously we use Lemma [7J] in- 
stead of Lemma 12.21 □ 

We now show an estimate for 6 defined in (|7.6[) . including exceptionally the 
exponential decay of the bound. 

Lemma 7.4. For the error 5 defined by (|7.6j) . we have 

\\6{t) || +h\\W5(t) || < Ch 2 e~ ct \\v h ||, fort>0, v h £ S h , with c> 0. 

Proof. Using the fact that Eh{t)Ah = —D t Eh(t), Lemmas 17.21 and 17.31 and the 
smoothing property (|2.2I) . we find this time taking into account the exponential 
decay of Eh(t) and Uh(t) for large t, 

\\S(t)\\ + h\\VS(t)\\ < J (\\E' h (t - s)Q h u h {s)\\ + h\\VE h (t - S )A h QhU h ( S )\\)ds 

<cj^ (t-sy^e-^-^^WQhUhis^l + hllAhQhUhis^ds 

< Ch 2 [\t-s)- x l 2 e- c ^\\Vu h {s)\\ds 
Jo 

<Ch 2 f {t-s)- 1 ' 2 e-< t - s h- 1 ' 2 e- cs ds\\v h \\=Ch 2 e- ct \\v h \\, 
Jo 

which is the desired result. □ 
We are now ready for the error estimates for the solution of (|7.2p . 
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Theorem 7.1. Let u and Uh be the solutions of ()1.19|) and (|7.2j) . Then for t > 0, 

H~ M /.si, ^ /C/J 2 |w| 2 , 7/||u h -u|| < C/i 2 |w| 2 , 

\Ch z t 7 |v|i, ifv h = P h v and \y7PhV\\ < C\v\i. 

Further, the estimates for the gradient of the error of Theorem \3.5\ remain valid. 

Proof. As in Section [31 it is suffices to estimate 8 = Uh — Uh- Using the splitting 
(|7.6p , 8 = 5 + 8, the term 8 is easily bounded by Lemma T7.41 and 8 is bounded as 
in Theorems 13.11 and 13.21 now applying Lemmas 17.21 and 17.31 □ 

Turning to nonsmooth initial data, we begin with the following lemma. 

Lemma 7.5. Let u and Uh be the solutions of (11.19P and (17. 2p . Then for t > 

\\u h (t) - u(t) - E h {t)A h Q h v h \\ < Ch 2 t- l \\v\\, ifv h = P h v. 

Proof. Using Lemma 17.41 for 8, it remains to bound S(t) — Eh(t)AhQhVh, which as 
for Lemma |3~TI is done as in [5J Theorem 4.1]. □ 

The following is now our nonsmooth data error estimate. Its proof is an obvious 
modification of that of Theorem 13.31 using Lemmas 17.21 17.31 and 17.51 

Theorem 7.2. Let u and Uh be the solutions of (|1.19|) and (|7.2j) . and let Qh be 
defined by (|7.5|) . Then, if (|1.17p holds, we have 

||5/,(i)-«WH^ a * _1 IMI. ifv h =P h v, fort>0. 

Condition 11.1 7(1 on Qh is again satisfied for symmetric meshes: 

Theorem 7.3. For {Th} symmetric, (|1.17jl holds for Qh defined by (|7.5p . 



Proof. We follow the steps in the proof of Theorem l4.ll For given x S Sh we define 
ip = ip x G H 1 as the solution of the Dirichlet problem A(p = \ in f — on dil. 
Since f2 is convex we have tp € H 2 and \1p\2 < C||xll- F° r ^ G ^/u we have 

llr) ,n {Qhip,x) a{Q h ^,p) 
\\Qhi>\\ = sup — — — = sup — 

xeSh llxll x es h \\x\\ 

\a(Q h ip,(p - I h ip)\ \a(Q h ip,I h <p)\ 
< sup — h sup — =1 + LI. 

xes h llxll x es h llxll 

By the obvious error estimate for Lh and Lemma 17.31 with p = 0, we get 

\!\<Ch sup m f}M 2 <Ch*W\\. 
xes h llxll 

To estimate II, we rewrite the numerator in the form 

a(Q h 'i/},Ih<p) = -s h (Qhip,Lhp) +a h (Qhi',JhLh<p) = ih + Hi- 
Id. order to complete the proof it suffices to show that 

+ ii 2 \ < Ch 2 llxll || VII- 
Using Lemmas 17.11 and 17.31 we obtain 

Kill < Ch 2 \\VQhi>\\ IIV/^II < C^UVQ^II \\<p\\ m < Ch 2 U\\ \\ x \\- 
Also, employing (|7.5p and (14.ip we get 

"2 = e h (ip,I h (p) = [ip,M h I h ip]. 
Since the family {Th} is symmetric, (|4. 12[) shows the required bound for ii 2 . □ 



ERROR ESTIMATES FOR THE FVEM 



27 



The results of Theorems 14.21 and 14.31 for our less restrictive assumptions on the 
family {Th} also remain valid, with the obvious modified proofs. 

The above results for the spatially semidiscrete finite volume method (|1.26p 
extend in the obvious way to the fully discrete backward Euler method (|6.ip and 
the Crank-Nicolson method (|6.3I) , with —A/, replaced by Ah, so that Theorems 
16. 1116.51 remain literally valid in the general case. 
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